Integrability and action operators in quantum Hamiltonian systems 
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For a (classically) integrable quantum mechanical system with two degrees of freedom, the func- 
tional dependence H = Hq{Ji, J2) of the Hamiltonian operator on the action operators is analyzed 
and compared with the corresponding functional relationship H{pi,qi;p2,q2) = Hc{Ji, J2) in the 
classical limit of that system. The former is shown to converge toward the latter in some asymp- 
totic regime associated with the classical limit, but the convergence is, in general, non-uniform. The 
existence of the function H = Hq{Ji, J2) in the integrable regime of a parametric quantum system 
explains empirical results for the dimensionality of manifolds in parameter space on which at least 
two levels are degenerate. The comparative analysis is carried out for an integrable one-parameter 
two-spin model. Additional results presented for the (integrable) circular billiard model illuminate 
the same conclusions from a different angle. 



I. INTRODUCTION 

A conspicuous phenomenological discriminant be- 
tween quantized integrable and nonintegrable parametric 
Hamiltonian systems with two or more degrees of free- 
dom is the occurrence or prohibition of level crossings be- 
tween states within the same invariant Hilbert subspace 
of the underlying symmetry groupuu'u Consider a quan- 
tum system with d continuous parameters whose classical 
counterpart is integrable on a manifold of dimensionality 
dj < d in parameter space. Empirical evidence suggests 
that level crossings occur on {dj — l)-dimensional man- 
ifolds which arCrembedded in the integrability manifold. 
A recent study,El which investigated this issue system- 
atically, showed for a two-spin model with d = 6 and 
di — 5 that the level crossing manifolds are, in fact, 
four-dimensional and that they are all confined to the 
five-dimensional integrability manifold. It showed, more- 
over, that the (classical) integrability manifold can be 
reconstructed from the (intrinsically quantum mechani- 
cal) level crossing manifolds. 

The focus of the present study is to illuminate the 
natural cause underlying this characteristic relationship 
between level crossing manifolds and integrability man- 
ifolds. We attribute this relationship to the presence of 
action operators as constituent elements of the Hamilto- 
nian operator for integrable quantum systems. 

The textbook solution of an integrable classical dy- 
namical system with two degrees of freedom, specified by 
an analytic function H{pi, gi;p2, 92) of canonical coordi- 
nates, is to transform the Hamiltonian into a function of 
two action coordinates: H = Hc{Ji, J2)- The canoni- 
cal transformation (j>i,qi) {Ji,Oi),i = 1,2 to action- 
angle coordinates amounts to a solution of the dynamical 
problem because it transforms Hamilton's equations of 
motion, pi — —dH/dqt, qi = dH/dpi, generically a set 
of coupled nonlinear differential equations, into Ji = 0, 
9i — dHc/dJi = LOi with the solutions Ji = const, 

e,{t) ^ uj,t + ef\ 

This solution is guaranteed whenever a second integral 
of the motion can be found, i.e. an analytic function 



/(pi, (7i;p2, 92) which is functionally independent of H 
and has a vanishing Poisson bracket with H: dl/dt — 
{iJ, /} = 0. Deriving the expressions i?c(Ji, J2) and 
/c(Ji, J2) from H and / requires the use of separable 
canonical coordinates. Finding separable coordinates can 
be a difficult task even if the second invariant is known. 

The functions Hc{Ji,J2) and Ic{Ji,J2) establish a 
pivotal link between an integrable classical system and 
a quantized version of it. Semiclassical quantization de- 
rives its raison d'etre from the obvious fact that quantiz- 
ing a functional relation is much less problematic if it in- 
volves only quantities such as H, /, Ji, J2 whose quantum 
counterparts are guaranteed to be commuting operators. 



II. QUANTUM VERSUS QUANTIZED 

In the context of this study, it is useful to distinguish 
and compare three versions of the same model system: (i) 
the quantum version, (ii) the classical version, and (iii) 
the (semiclassically) quantized version. 

The primary version is the quantum model, specified 
by the Hamiltonian expressed as an operator valued func- 
tion of a set of dynamical variables (position, momentum, 
spin, . . . ) The commutation relations of these operators 
and the metric of the associated Hilbert space along with 
the rules of quantum mechanics then determine, via the 
Heisenberg equation of motion, the time evolution of any 
observable quantity of interest. 

The classical limit converts the Hamiltonian operator 
into the classical energy function, the commutator alge- 
bra of dynamical variables into the symplectic structure 
(the fundamental Poisson brackets), and the Heisenberg 
equation of motion for any operator into the Hamilton 
equation of motion for the corresponding classical quan- 
tity. These quantities, in turn, enable us to express the 
energy function as a classical Hamiltonian, i.e. as a func- 
tion of canonical coordinates. 

The quantization of a classical Hamiltonian system re- 
quires a prescription for translating the functional re- 
lations between classical dynamical variables into func- 
tional relations between corresponding operators. Semi- 



2 



classical quantization is one neat and clean procedure 
applicable to all integrable classical systems. It bor- 
rows from classical mechanics the functional dependence, 
H = Hc{Ji, J2), of the Hamiltonian on the action oper- 
ators and postulates that the eigenvalue spectrum of the 
latter consists of equidistant levels spaced by h: 



1 



1,2 



(1) 



with integer n^. The (integer) Maslov indices ai are 
determined by the topology of the classical trajectories 
in phase space.Ll Semiclassical quantization thus makes 
specific predictions for the energy level spectrum of the 
quantized version of the model system at hand. 

In general, the (semiclassically) quantized and the (pri- 
mary) quantum energy level spectra of one and the same 
integrable model system do not coincide. The relation- 
ship between the two spectra will be investigated in 
Sec. [II for an integrable two-spin model and in Sec. |^ 



for the (integrable) circular billiard model. 



III. TWO-SPIN MODEL 

We consider two quantum spins Si, S2 of equal length 
\/o{p -I- 1) (cr = i, 1, |, . . . ) interacting via a uniaxially 
symmetric exchange interaction^ 



(2) 



The second integral of the motion, which follows from 
Noether's theorem, is 



(3) 



z(t) ^ 



-s(cos79i — COS 1^2) = zosn(p<, zo/a), 



C(t) = tan(^i - (^2) 
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where sn(p, x), cn(p, a;), dn(p, x) are Jacobian .elliptic 
functions and K(p) is a complete elliptic integral.cl 

For the case k = 1 with higher rotational symmetry, 
considerable simplifications occur in the classical time 
evolution. Both spins precess uniformly about the di- 
rection of the conserved vector Sj^ = Si S2, and the 
precession rate is w — \St\ for both spins. Equations (^) 
for the classical actions become 



Ji 
J2 



2M^, 
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z{t) — zosinai, C,{t) = 



(l + C2)(S + M|-z2) 
azQ cos at 



(5a) 
,(5b) 
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E), 



and can be evaluated in closed form: 
Ji = 2M„ 



(6a) 

-^2{s-^-E) + {s ~ M,)sgn{s^ -E- 2sAQ 



[s + M^)sgn{s^ -E + 2sM^). 



(6b) 



Inverting relations (||) yields a degree-two polynomial de- 
pendence of E, Mz on Ji , J2 : 



In the classical limit /i — s- 0, ti — > 00, h^J a{a + 1) = s, 
the operators S.^ turn into 3-component vectors. Si = 
s(sinz9i cos (/9i, sin i?; sin c^;, cosiJi), and Eq. (|^) then de- 
scribes the energy function of an autonomous Hamilto- 
nian system with two degrees of freedom aBfi canonical 
coordinates pi = scosdi, qi — cpi, i = 1, 2.1 



/c(Ji, J2) = M, - ^Ji, 

where Ic = J2 ^ \ Ji \ if s| Ji| > — E and Ic 
s\Ji\ <s^-E. 



(7a) 
(7b) 
2s ~ J2 if 



A. Classical actions 

Generically, the classical time evolution of this system 
is nonlinear and quasiperiodic. In the parameter range 
< K < 1, the following relation between the integrals 
of the motion H = E (energy), / = Mz (magnetization) 
and a set of classical actions Ji , J2 can be inferred from 
the exact solutioniEl 



B. Quantum actions 

For the case k = 1, the exact quantum spectrum fol- 
lows directly from the higher rotational symmetry of H: 

{H)q = h^a{a + 1) - + 1), {Mz)q = Jm, (8) 

where Z = 0, 1, . . . , 2(7 is the quantum number of the total 
spin and m = —I, —l + l,... ,+l that of its z-component. 
One set of quantum actions (|l]) has eigenvaluesll2l 

(J,)/ri= = -a,-c7-M,... ,-t-a, (9) 



3 



which are related to l,m as follows: 



Q 



r 

tQ 



a -I, J^=cr-l-m (m < 0), (10a) 



(T-l + m, J§^a~l (m>0). (10b) 



The two quantum invariants expressed as explicit func- 
tions of action operators then read 

Hq{JiJ2)^H = iftV(f7 + l) + imin(Ji,J2) 

X [ft(2cr + l) -min(Ji,J2)], (11a) 



/q(Ji, J2) = Af, = ^iJi~J2), 



(lib) 



where min(Ji, J2) selects the action operator with the 
smaller eigenvalue. 

While the functional dependence in ( pi] ) is again de- 
scribed by a degree-two polynomial, it is different from 
the functional dependence ^ found classically. The for- 
mer cannot be reconciled with the latter by any canonical 
transformation, nor does the quantum spectrum converge 
uniformly toward the classical spectrum for cr —> 00, as 



we shall see in Sec. Ill C 1 . 

For the cases < k < 1 we must calculate the 
(2a + 1)^ eigenvalues of the two quantum invariants 
H, Mz by numerical diagonalization of H in the 4cr + 1 
invariant subspaces of M^. From the numerical data 
for (H), {Mz), we can infer the correct assignment 
of action quantum numbers (Ji) /h to eigenstates by 
smoothly connecting the spectrum in parameter space to 
the known relations (^1|) for k = 1. The resulting data 
for Hq{Ji, J2), /q(Ji, J2) can then be compared with the 
(semiclassically quantized) inverse classical relations (jj), 
Hc(,Ji, J2), Iq(,Ji, J2), to high precision albeit not ana- 
lytically as in the case k = 1. Numerical results will be 



presented in Sec. IIIC2 



be compared to the semiclassical expressions Hc{ji, ^2), 
/c(Ji, 0/2) inferred from the classical relations (^ with 
quantum actions It turns out to be more practi- 

cal to perform the comparison for the inverse functional 
relations. We substitute a{(j + 1) for and the exact 
eigenvalues (||) for E, Mz into the classical expressions 
(^ . The result is a set of non-integer valued semiclassi- 
cal action quantum numbers 
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(12a) 

TO / = 

2y^a{a + l)- y /l{l + l) \m\ < too (12b) 
|to| - y^l{l + 1) |m| > Too, 



where Too = l{l+l)/2y^a{a + 1). An optimal match with 
the quantum actions (|l^) can be achieved if we subject 
(O) to two successive canonical transformations: 



T C 



C 



2^a{a + l) 



■^2 



J? 
J? 



Jl 



J2 



J2^ - 2^/^(0- + l) + a+i 



< 
>0, 



Ji^' < 



T c" 

J2 

T C 
J2 



2^a(a + l) + a + Ji^' + i Ji^' > 



2^cr(fT + l) + cr- 
J2<^'- 2^/^(0- + l) + cr- 



Ji^ < 
Ji^'>0. 



We thus arrive at the expressions 



J/ 



J2^ - 



m = l = 

TO < (13a) 



cr - ^/(/ + 1) + i + TO TO > 0, 



1^0 



a - y /l{l + l) -m+^ TO<0 (13b) 
cr - + 1) + i TO > 0. 



C. Quantum corrections to quantized actions 

In some simple applications, the functions Hq,Iq are 
identical to the functions Hc,Ic- Hence there are no 
such quantum corrections. If we take, for example, the 
two-spin model H = — S'fS'f, then both classical invari- 
ants E, Mz depend solely on the canonical momenta, and 
the latter are identified to be actions: pi = Ji. Hence we 
have E — — J1J2, A/^ ~ \{Ji + J2)i which, upon semi- 
classical quantization with {Ji) /h — —a, — cr-|- 1, . . . , -|-cr, 
yields the exact quantum eigenvalue spectrum. This sit- 
uation is exceptional. For all cases of (g) with < k < 1, 
quantum corrections do exist. 



1. Exact results for k = 1 

For the parameter setting k = 1, the functions 
Hq{Ji, J2), Iq{Ji,J2) as given by expressions (111]) are to 



The deviations of the non-integer valued Jl , J2 from 
the integer valued J^ , J^ then describe the quantum cor- 
rections to the semiclassical actions. 

Using + 1) — i = ? + 0{l~^), we see at once that 
the genuinely quantum mechanical relations ( |l0| ) and the 
semiclassical relations (O) are asymptotically equivalent 
at low energies (large f) for ct ^ cxi. At high energies 
(small I), on the other hand, the two relations remain 
distinct no matter how large we choose the value of the 
spin quantum number cr. 

To set the stage for the cases < k < 1, we plot in 
Figs, ^a) and ||(a) the eigenvalues of H versus those of 
Mz in representations with spin quantum numbers cr = 2 
and cr = 4, respectively. The patterns of regularity and 
similarity in the arrays of points are a direct consequence 
of the smooth functional relations Hq {Ji, J2), Iq{Ji, J2)- 
The map {{H), {Mz)) {jf, J^) from the plane of in- 
variants to the action plane is provided by Eqs. ( p^ ) 
and produces the triangles in Figs. |](b) andH(b). These 
points form a perfect lattice with unit spacing. 
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FIG. 1: (a) Eigenvalue (H) (energy) versus eigenvalue (M^) 
(magnetization) as given in Eqs. (^) of all eigenstates of 
Hamiltonian (0) with k = 1, a = 2. (b) The full triangles are 
the quantum images (J^, J^) of these eigenstates in the ac- 
tion plane as provided by Eqs. (jlo|) . The open circles are the 
semiclassical images (Ji'^ 
with = ct((7 + 1). 





J1 



J2 ) as provided by Eqs. 



If we use instead the map ( p^ ) provided by semi- 
classical quantization, we obtain the array of open cir- 
cles in Fig. |l](b) and Fig. ||(b). The bonds shown in 
parts (a) and (b) of the two graphs correspond to each 
other. The distortion in the lattice of circles relative 
to the perfect lattice of triangles is a graphical repre- 
sentation of the quantum corrections in the functions 
Hq{Ji, J2), /q(Ji, J2) relative to the semiclassical func- 
tions Hc{Ji, J2), Ic{JiiJ2)- It visually confirms what 
we have already concluded from comparing ( p^ ) and (13), 
namely that the deviations die out at low energies (lower 
left area) but persist at high energies (upper right area) 
for a ^ 00. A useful measure of the leading quantum 
correction to the semiclassical relation Hc{Ji, J2) is the 
quantity ctA J, where 



A J = v/(AJi)2 + (AJ2)2, A J, = Jf^-J^ 



(14) 



represents the distance between the triangles and circles 
on corresponding array sites in Figs. |](b) and 0(b). From 



FIG. 2: Plot of the same quantities as in Fig. |l| but for spin 
quantum number a = 4. 



Eqs. (nq) and (Hi) we obtain 



AJ = 



I/V2 1^0 



(15) 



The dependence of aAJ on J{^, thus represents the 
1/(7 quantum correction to the semiclassically quantized 
actions. It has an inverse first power divergence in one 
corner of the action plane for energy levels at the upper 
threshold of the spectrum: a A J ^ [4^/2(1 / a)]~^ . For 
states with l/a <^ 1 the leading quantum correction is of 
0(1). In this part of the spectrum, semiclassical quanti- 
zation remains inadequate no matter how large we choose 
the spin quantum number a. 

The state with the largest quantum correction to semi- 
classical quantization is the singlet combination of the 
two spins. This state or any nearby state in the action 
plane have no proper semiclassical representation. 



2. Numerical results for < k < 1 

Here we use the same graphical representation even 
though we must rely on the results of a numerical di- 
agonalization for the energy eigenvalues. At k < 1 we 
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FIG. 3: (a) Eigenvalue (H) (energy) versus eigenvalue (M^) 
(magnetization) of the {2a -f 1)'^ = 25 eigenstates of the two- 
spin model (^ with k — 0.1 for a — 2. Data from a numer- 
ical diagonalization. (b) The full triangles are the eigenval- 
ues j'^ = {Ji)/h of the action operators, the images of the 
inverted functions Hq{Ji, J2), Iq{Ji, Ji)- The open circles 
are the semiclassical images (Ji^",J2^") fromEqs. (g) with 
— a{a+l), the images of the inverted functions Hc{Ji, J2), 



observe that certain features of the quantum invariants 
change quahtatively because the rotational symmetry of 
H has been reduced, whereas other features remain qual- 
itatively the same because the integrability of the model 
has not been destroyed. 

In Figs. ||(a) and ^a) we have plotted the eigenvalues 
(if), (Mz) of the two quantum invariants versus each 
other at K = 0.1 for (T = 2 and = 4, respectively. 
Again the data points display regular patterns. They 
evolve from the patterns shown in Figs. [l](a) and ^(a) by 
smooth deformation of the lines of bonds as the value of 
K is lowered gradually. The lower symmetry removes the 
level degeneracies pertaining to the strings of horizontal 
bonds in Figs. |l|(a) and ||(a). 

When we substitute the eigenvalues (H) and {Mz) 
from the numerical diagonalization into the exact ex- 
pression (||) for the classical actions and subject the re- 
sulting set of discrete values Jf to the transformations 
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FIG. 4: Plot of the same quantities as in Fig. |^ but for spin 
quantum number a = 4. 



jf — > Ji^' Ji^" , we obtain arrays of points in the 
form of distorted lattices as illustrated by the open circles 
in Figs, ^(b) and ^(b) for the two examples at hand. The 
deviations of these data points from the sites of a per- 
fect lattice (marked by triangles) then again represent 
the quantum corrections to the (semiclassically) quan- 
tized actions. The patterns in Figs, ^(b) and |j(b) are 
also connected to those in Figs. |l|(b) and g(b) by smooth 
deformation of the lines of bonds upon gradual variation 
of the parameter n. 

A closer look at the 1/ a quantum correction is afforded 
if we plot the scaled distance crA J versus the scaled ac- 
tion quantum numbers j a and j a for a system with 
many more levels (cr = 40) . A contour plot of the result- 
ing landscape is shown in Fig. |5|. Convergence of ctAJ 
toward a smooth function of Ji/o, j'^ /a is almost uni- 
form. In the case k = 0.1 considered here, there are two 
points (as opposed to a single corner point at k = 1), 
where the 1/cr correction diverges. The data points crAJ 
closest to these locations again tend to grow ex a. 

The two sharply peaked maxima in the landscape of 
Fig. ^ will merge into a single divergence as ct — *■ cxd. 
At this point in the action plane, the leading quantum 
correction to semiclassical quantization is again of 0(1). 
Its location in the action plane does, however, no longer 



-1 -0.5 0.5 1 

J?/a 



-1. -0.5 0.5 1. 

J°/o 



FIG. 5: Scaled distance crAJ for o — 40, k = 0.1 between the 
images of the inverted functions Hq{Ji, J2), Iq{Ji,J2) and 
the images of the inverted functions Hc{Ji, J 2), Ic{Ji, J2)- 



FIG. 6: Scaled distance crA J for a = 40, k = 0.5 between the 
images of the inverted functions Hq{J\, J2), /q(Ji, and 
the images of the inverted functions Hc{Ji, J 2), Ic{Ji, J2)- 



coincide with an extremum in the energy level spectrum. 
The divergence in crAJ occurs at energy E = ks^ (for 
CT — > 00), where the classical equations of motion have a 
fixed point. For eigenstates with action quantum num- 
bers in the vicinity of this point, quantum effects persist 
no matter how large a is made. 

One point in the action plane where ctA J diverges ex- 
ists throughout the regime < k < I. With k increasing 
from zero, the singularity moves gradually toward one 
corner of the action plane, and the energy of the state 
pertaining to those action coordinates moves toward the 
upper threshold of the spectrum. This trend is indicated 
in Fig. |[ which shows the l/cr-landscape for k ~ 0.5. 
The endpoint of this gr adual shift, the case k — 1, was 
described in Sec. |III C l| . 

The asymptotic landscape for cr ^ 00 to which the 
graphs in Figs. ^ and |6| converge almost everywhere can 
now be used as the reference frame for the higher-order 
quantum corrections. The deviations of the data points 
from this new reference, appropriately scaled, will pro- 
duce another landscape, representing theJ/cr^ correction 
to the semiclassically quantized actions.llil 

We consider the line J? = Ji — <j/2 for this purpose. 
In the main plot of Fig. ^ we show the 1/a corrections 
crAJ along this line for cr = 4, 8, 16, 32. Also shown are 
data for a — 1600, which are very close to the asymp- 
totic values for the 1/a correction and now serve as the 
reference line for the l/a^ corrections. 

In the inset to Fig. |^ we have plotted the scaled devi- 
ations of the a = A,8, 16, 32 data from the new reference 
line. The results suggest that these data again converge 



toward a line, which will then be the reference line for 
1 /a^ corrections. Like the reference line in the main plot 
of panel (a) [panel (b)], which is embedded in the land- 
scape Fig. 1^ [Fig. ^, the new reference line will be em- 
bedded in a landscape representing the l/cr^ quantum 
corrections to semiclassical quantization over the entire 
action plane. 

The point to be emphasized here is not the exact 
shape of the landscapes that represent successive orders 
of quantum corrections to the semiclassically quantized 
actions, but that such corrections exist and that the lead- 
ing term may be of 0(1) at special points rather than of 
0{a~^) as might be expected. 

IV. CIRCULAR BILLIARD 

In the second application we consider a particle of mass 
m that is free to move two-dimensionally across a circular 
area of radius R. The classical Hamiltonian expressed in 
polar canonical coordinates reads 

H{pr,r;p^,^9)^^ + ^ + V{r), (16) 

where V{r) is a hard- wall potential that confines the par- 
ticle to r < R. |— I 

In a recent study, Ree and Reichlli3 analyzed this sys- 
tem classically and quantum mechanically as an inte- 
grable limiting case of the circular billiard with a straight 
cut. In general, the cut renders the classical time evo- 
lution chaotic. Here we use some results of Ref. O to 
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FIG. 8: Eigenvalue (H) (energy) versus eigenvalue (L) (an- 
gular momentum) as given in Eq. ( p^ of the eigenstates near 
the bottom of the spectrum of the circular billiard model. 



with = 2mE/h^ and Dirichlct boundary conditions 
are known. The exact expressions for the two quantum 
invariants H (energy) and L (angular momentum) are 



FIG. 7: Dependence of the scaled distance aAJ for (a) 
K — 0.1, (b) K = 0.5 between the images of the inverted 
functions Hq{Ji, J2), Iq[J\,J2) and the images of the in- 
verted functions Hc{Ji, J2), Ic{Ji, Ji)- Shown are data for 
cr = 4 (squares), a = 8 (circles), a = IQ (triangles), cr = 32 
(pentagons), and a — 1600 (solid line). Inset: Scaled devia- 
tion alaAr^f - aAJ"] of the a = 4,6,8, 16 data from the 
reference line (ct = 1600 data). 



investigate the functional dependence of the circular bil- 
liard Hamiltonian on the actions quantum mechanically 
and semiclassically for comparison with the two-spin re- 
sults presented previously. 

Integrability of the circular billiard model is guaran- 
teed by the conservation of angular momentum L = p^. 
The canonical transformation to action-angle coordinates 
produces the following relations between the integrals of 
the motion E, L and the two-action variables: 



(H) = 



"Ik 



Ji = L, 



y/2mE 



(17a) 

:arccos (^j^ i (17b) 



where x = L"^ /2mE. The eigenfunctions of the circular 
billiard, i.e. the solutions of 



r dr 



1 

J.2 Q^2 



(18) 



{L) = ±in, 



(19) 



where I — 0,1,2, . . . and aik is the fc*^ zero (fc = 1,2,...) 
of the Bessel function J; (a;). 

One major distinction between the circular billiard 
model and the two-spin model is that all invariant Hilbert 
subspaces are infinite-dimensional in the former and 
finite-dimensional in the latter. The energy has no upper 
bound in the circular billiard and the angular momentum 
neither upper nor lower bound. 

In Fig. 9 we have plotted the eigenvalues (H) ver- 
sus (L) of the two quantum invariants near the bottom 
of the level spectrum. As in the two-spin model, the 
regular pattern of points is a signature of quantum in- 
tegrability. In both models the points tend to become 
displaced irregulajly when nonintegrable perturbations 
are introduced.EHl3 

The integers k,l in (|l9|) can be identified as the eigen- 
values (in units of h) of a set of quantum actions: 



(Ji) = hi, {.h) = h{k - i). 



(20) 



The shift in the second expression is dictated by a Maslov 
index ai = 1 (see Sec. II) The results ( p^ ) com- 
bined with ( pO| ) thus define specific functional relations 
Hq{Ji, J2), Iq{Ji,J2) between quantum invariants and 
quantum actions. They are to be compared with the 
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functional relations Hc{Ji, J2), Ic{Ji, 'h) as defined by 
( |l7| ) combined with (|20|). 

For a graphical representation of the quantum cor- 
rections to semiclassical quantization, we proceed as in 
Sec. In Fig. ^ we plot AJ2 ee | J;^ - | versus k and 



^ and Jo is the value of (17b) when 



Z, where J2 = k — j anu ^2 

the exact eigenvalues ( p^ ) for the quantum invariants are 
substituted into the expression. 

We observe a landscape in the form of a sloped ridge 
centered at Z = 0. The largest quantum correction to 
semiclassical quantization pertains to the ground state 
(with k — 1,1 — 0). The plot suggests that the quan- 
tum corrections die out for large k. This is confirmed by 
substitution of the asymptotic expression for k ^ Ip 

4/2 7 1 

c^ik-^p-^ + Oir'), f3 = k+---, (21) 



into (|l|) for use in (|l7b| ): 
1 1 



J2(z,fc) ~ h 



fc - ^ + 0(fc-2) 

4 87r2fc ^ ' 



, k->l. (22) 



The quantum corrections also decrease with increasing 1 1 \ 
at fixed k, but not all the way to zero. To demonstrate 
this fqs fc = 1, we use the asymptotic expression for I ^ 
k = lE 



an ~ |^|+Ci|/|i/3 + C2|r' 



/3 



(23) 



with Ci ~ 1.8558, C2 ~ 1.033 for use in (|T|). When 
substituted into (|l7bD we obtain the asymptotic value 



.h{l, 1) = (n/37r)(2C7i)3/2 -I- 0(1/1-2/3), (24) 

which deviates from the reference value h{l — by 
roughly one percent. The conclusion is that the semi- 
classical regime of the circular billiard is restricted to 
states with k ^ I. It does not include, for example, an 
states along the lowest branch {k — 1) shown in Fig 
no matter how large the energy of the state becomes with 
increasing |/|. 

V. CONCLUSION 

In this study we have investigated a key signature of 
quantum integrability in systems with two degrees of free- 
dom, namely the functional dependence of the Hamilto- 
nian H and the second integral of the motion / on two 
action operators Ji , J2 . 



The results presented in Sees. ID and IV for the (semi- 
classically) quantized and the (primary) quantum en- 
ergy level spectra of two integrable model systems sug- 
gest the following interpretation, which is consistent with 
the conclusicps inferred from an entirely different line 
of reasoningil3 (i) Quantum integrability implies that 
the Hamiltonian can be expressed as an operator val- 
ued function of the actions: H = Hq{Ji, J2), where 




FIG. 9: Quantum corrections to the semiclassical prediction 
for the energy eigenvalues of the circular billiard model. Plot- 
ted is the deviation AJ2 = \J? — JI^I, where J? = k — 



Jk = 'h/h as determined by (17b) with E = {H),L = {L) 



J and 



substituted from ( p^ 



the eigenvalue spectrum of the action operators is of the 
form dl^) . (ii) This function is different from the function 
Hc{Jit J2) inferred via semiclassical quantization from 
the solution of the classical dynamical problem, (iii) 
In some asymptotic regime associated with the classi- 
cal limit the function Hq{Ji, J2) converges, if properly 
scaled, toward the function Hc{Ji, J2), but the conver- 
gence need not be uniform, (iv) For the second integral 
of the motion, which (classically) guarantees integrabil- 
ity, there exist functions /q(Ji, J2) and Ic{JiiJ2) with 
analogous properties. 

The existence of action operators as constituent ele- 
ments of all quantum invariants in integrable model sys- 
tems is a key property necessary to explain the dimen- 
sionality of level crossing manifolds relative to the di- 
mensionality of integrability manifolds in the parameter 
space of model systems with parametric integrability con- 
ditions. On the ^/-dimensional integrability manifold in 
the parameter space of a given model system, both func- 
tions Hq{Ji, J2) and /q(Ji, J2) will then depend con- 
tinuously on these parameters. The quantum eigenvalue 
spectrum on the integrability manifold is determined by 
{H)q — Hq{{Ji), ( J2)) and can be interpreted as a set of 
continuous functions of the Hamiltonian parameters sub- 
ject to the constraints imposed by the integrability condi- 
tion. The level crossings, which occur at the intersections 
of the graphs of any two members from the set of func- 
tions are then naturally confined to (c?/ — l)-dimensional 
manifolds and are naturally embedded in the integmbil- 
ity manifold, in agreement_with empirical evidence .□ 

In a companion paper ,0 we have investigated how 
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the existence of Hq{Ji, J2), lQ{ji, within the five- 
dimensional integrabihty manifold of a six-parameter 
two-spin model affects the properties of quantum invari- 
ants and what impact on the same quantities the nonex- 
istence of Hq{Ji, J2), Iq{Ji, J2) elsewhere in parameter 
space has. 
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